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ABSTRACT 



This paper is an application of the Piecewise-Sinusoidal 
Reaction Matching Technique (PSRMT) to linear, parallel 
thin-wire dipoles. The "method of moments" if briefly 
outlined and the fundamental electromagnetic principles 
of reciprocity, surface equivalence and reaction are dis- 
cussed. Through a logical combination of these concepts, 
PSRMT is developed and applied to thin-wire dipoles radiating 
in free space due to a delta function voltage generator at 
their center. The driving point impedance and the current 
distribution are calculated and compared to results obtained 
from other independent theoretical techniques. 

The major contribution of this investigation is not that 
of presenting a new technique for calculating the input im- 
pedance and current distributions on dipole antennas. The 
study is primarily concerned with comparing other well-known 
techniques with that of PSRMT. Emphasis is placed on the 
ease and convenience of expressing the mutual impedance 
integrals of coplanar dipoles in closed form. 
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I . INTRODUCTION 



A. DESCRIPTION OP ANTENNA PROBLEM UNDER INVESTIGATION 

This paper is an investigation of the efficiency of 
the Piecewise-Sinusoidal Reaction Matching Technique (PSRMT). 
In the solution of problems using the "method of moments" 

[3], computation time on the computer depends primarily 
upon the amount of computational effort required to generate 
the elements in the generalized impedance matrix. Usually 
this relates to the amount of numerical integration required 
to calculate the Z^'s. However, by employing piecewise- 
sinusoidal basis functions to obtain the impedance matrix, 
closed form expressions can be obtained for the integrations. 
Also, the size of the impedance matrix required to adequately 
solve the antenna problem is relatively small when compared 
to the use of other basis and weighting functions. For 
these reasons PSRMT will be applied and the results compared 
to other techniques. 

The antenna configuration under investigation will be 
restricted to a linear, thin-wire dipole. The antenna is 
assumed to be radiating in free space and is driven at the 
center by a delta function voltage generator. In this 
particular application of PSRMT, all segments of the fila- 
mentary dipoles are parallel and of equal length. This is 
not a necessary condition, but one of convenience. The 
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driving point impedance and the current distributions 
will be calculated by an application of PSRMT. 

B. BRIEF HISTORY OF RELATED WORK 

In the past much work has been done in the area of 
determining the mutual impedance between two elements with 
an assumed sinusoidal current distribution. Carter [5] 
has presented the mutual impedance between two wires in 
echelon where each antenna is an odd number of half-wave- 
lengths long. The mutual impedance expression between 
two identical nonstaggered parallel center-fed antennas 
of Brown [1] is of historical renown. The work of Brown 
was carried one step further by Cox [6] who developed the 
mutual impedance between parallel antennas of unequal lengths. 

In 1957, H. E. King published a paper [10] that developed 
an expression for the mutual impedance between two staggered 
parallel, center-fed, infinitely thin antennas of unequal 
lengths. Lewin [12] analyzed the mutual impedance of two 
coplanar half-wave dipoles, in arbitrary positions. Finally 
Lewin’ s work was generalized by Richmond and Geary [10] 
for any length antenna. 

All of the above studies were restricted to assumed 
sinusoidal dipole current distributions. Without assuming 
a current distribution upon an antenna the current distri- 
bution may be found by breaking the antenna up into several 
segments, calculating the mutual impedance between each 
pair of elements in turn and generating an impedance matrix 
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that characterizes the current distribution. This calcula- 



tion of the mutual impedance is done in general by the 
expression published by Richmond and Geary. The technique 
crudely described is that of reaction matching. 

C. SCOPE AND LIMITATIONS OF THE STUDY 

This paper develops the multisegment reaction technique 
from fundamental electromagnetic theory. Although the 
technique is applied only to linear thin-wire dipoles of 
perfect conductivity, the procedure may be extended to 
more general and intricate geometries. To cite one of 
the many examples, PSRMT has been successfully employed 
by Travieso-Diaz [20] to the analysis of antennas and 
scatterers with arbitrary surface shape. 
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II. OUTLINE OF THE METHOD OF MOMENTS 



The following discussion portrays the major steps 
employed in the method of moments. The object is to 
reduce the functional equations of field theory to matrix 
equations which can be readily solved. The name "method 
of moments" has been given to the mathematical procedure 
for obtaining the matrix equations [8], 

A. FORMULATION OF THE OPERATOR EQUATION 

Consider only equations of the inhomogenous type 
L(f) = g, where L is a given operator, g is the excitation 
or source (known function), and f is the field or response 
(unknown to be determined). 

One must determine from physical considerations that 
the problem in question is deterministic, that is, if the 
solution for f when L and g are given is unique. The method 
of moments gives a general procedure for treating field 
problems, but the details of solution vary widely with 
the particular problem. The analysis to follow is that 
which is germane to the antenna problem under investigation. 

B. POSTULATION OF THE BASIS FUNCTION 

Consider the deterministic equation 

L(f) = g , (1) 
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where L is a linear operator, g is known, and f is to be 
determined. Let f be expanded in a series of functions or 
a linear combination of a set of N vectors 



where the c 's are constants that are not known but are 
n 

to be determined in the solution process. The f ' s are 
chosen to be linearly independent and are called expansion 
functions or basis functions. For exact solutions. Equation 
(2) is usually an infinite summation and f ' s form a complete 
set of basis functions. For approximate solutions, Equation 
(2) is usually a finite summation. Substituting Equation 
(2) into Equation (1), and using the linearity property 
of L, 



C. POSTULATION OF THE WEIGHTING FUNCTION 

It is assumed at this point that a suitable inner 
product <Cf,g> has been determined for the problem. The 
inner product is a scalar defined to satisfy 





( 2 ) 



y c L(f ) = g 

n n 
n 



( 3 ) 



<f,g> = 



^^f + £g,h> = «(<f,h> + ^<g,h> 



<f*,f> >0, if f / 0 

= 0, if f = 0 , 
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where o' and ^ are scalars, and * denotes complex conjugate 
[8]. The appropriate inner product for this particular 
antenna problem will be developed in the section of this 
paper entitled "Reaction Concept". Nov; define a set of M 
linearly independent weighting or testing vectors {wj}, and 
take the inner product of Equation ( 3 ) with each . The 

result is 



Z c /w ,Lf y = <w ,g"> ( 

n n \ m* n ^ m ,& ' v 

m= 1, 2, 3 , ••• M, where the N c^ ' s are unknown and all 
other terms are known [3]. 

D. SOLUTION OF THE LINEAR ALGEBRAIC SYSTEM 

Equation ( 4 ) can be written in matrix form as 

|\in] [ C n] [^n] 5 ^ 



where 




r i 






c i 






c 2 




-| 




v 7 2 ,g 




and 


iL 










n 






c 3 








7 3 ,g 


• 




• 



A consistent equation results by requiring M and N to 

be equal. That is, the basis set f and the testing set 

w each has N elements. The coefficient matrix is deter- 
m 

mined by inversion: 



N • W 

provided that the matrix ^ 1 is nonsingular. The 

solution for f is given by Equation (2): 



S. 



in 



( 6 ) 




( 7 ) 



This solution may be exact or approximate, depending upon 

the choice of F and w and the extent of the summation. 

n m 

The particular choice w = f is known as Galerkin's method. 

m n 

E. COMMENTARY ON THE METHOD OP MOMENTS 

One of the most critical tasks in any particular problem 
is the choice of the basis and testing functions, f n and 
w . The f's should be linearly independent and chosen so 



15 



that some superposition (see Equation (2)) can approximate 

f reasonably well. The w^’s should also be linearly 

independent and chosen so that the products ^w n ,s> 

depend on relatively independent properties of g. Obviously 

additional considerations affecting the choice of f and 

° n 

w are a.) the accuracy of solution desired, b.) the ease 
m 

of evaluation of the matrix elements, c.) the size of the 
matrix that can be inverted, and d. ) the realization of 
a well-conditioned matrix fT 1 [ 8 ] . 

I mn J 

F. PIECEWISE-SINUSOIDAL BASIS SETS 

The physical problem may be represented by many 

different operator equations, and a suitable one must be 

chosen. There are also an infinite number of sets of 

expansion functions f and testing functions w that may 

n m 

be chosen. 

Intuitively, the more nearly the expansion set resembles 
the actual form of the current distribution, the fewer the 
terms required to -accurately represent the current. In 
addition, certain types of basis functions may lessen many 
time-consuming, numerical operations incorporated in 
approximations. The piecewise-sinusoidal basis set exhibits 
these qualities. The objective is to represent the current 
distribution, over individual subsections of a domain with 
expansion functions that possess greater curve-fitting 
potential than do flat pulses and yet lead to well-conditioned 
matrices in approximate solution procedures. The application 
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of piecewise-sinusoids in Galerkin’s method and the 
method of moments to electromagnetic field problems is 
known as the reaction concept [3,17]. 
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III. THEORETICAL MOTIVATION 



The following discussion will develop the basic concept 
utilized in the formulation of the reaction matching 
technique of J. H. Richmond. The reaction concept of 
Rumsey [17], the reciprocity principle [9] and the surface 
equivalence theorem [ 18 ] were presented by Richards [13] and 
Travieso-Diaz [20]. The concepts of reciprocity, reaction 
and surface equivalence form the theoretical motivation 
for PSRMT; therefore they are included again here. With 
a logical combination of these elements, one can develop 
a general formulation of electromagnetic boundary-value 
problems, in particular, reaction matching. 

A. RECIPROCITY THEOREMS 

Consider a pair of time-harmonic sources of the same 
frequency w radiating in a linear, isotropic, but not 
necessarily homogenous medium with permittivity L and 
permeability f * . One of the sources, designated as source 
1, has impressed electric and magnetic current densities 
(J^, M^) which radiate fields (E^, H^) in the presence of 
a second source, source 2, acting as a receiver. Similarly, 
the current densities (Jg, Mg) source 2 set up fields 
(Eg, Hg) in the presence of source 1 acting as a receiver. 

The fields of these sources at every point in the 
medium are related by the point reciprocity theorem of 
Lorentz [9] as follows: 
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v . (E 1 X H 2 - E 2 X H x ) = (J 1 • E 2 - J 2 • E 1 ) 

- (f^ • H 2 - M 2 • H 1 ) 

( 8 ) 



Now let the sources be of finite extent, and assume all 
inhomogeneities in the medium to have a finite distribution. 
Define non-overlapping surfaces and S 2 enclosing the 
source distributions (J-^, M^) and (J 2 , M 2 ) respectively; 
all inhomogeneities are contained in the volumes and V 2 
bounded by and S 2 . Therefore, the medium outside 
and S 2 is homogeneous. Under these conditions, the 
reciprocity theorem of Carson [ J l] establishes a relation 
between the source distributions (J^, M^) and (J 2> M 2 ) 
and the fields (E^, H^) and (E 2 , H 2 ) of Equation (8) as 





• E ! " m 2 * ^Jdv 
(9) 



B. SURFACE EQUIVALENCE PRINCIPLE 

The surface equivalence principle can be briefly 
described as follows [13]. First, an original source 
distribution (J\, PL) set up fields (E, H) everywhere in 
some defined medium. Next, the original source is turned 
off, but simultaneously surface current densities (J , M ) 

o o 

are set up on a surface S which encloses the original 
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source and any inhomogeneities. If the surface current 



densities ( J , M ) are related to the original fields 
s s 



(E, H) by 



J g = n X H 



(10A) 



M s = E X n 



y 



(10B) 



A 



where n is the unit normal at S directed outward from S, 
then these surface current distributions will radiate the 
original fields (E, H) in the region exterior to the surface 
S and will generate a null field inside S. 

C. REACTION CONCEPT 

In 195*1 V. H. Rumsey published a paper [17] in which 
he introduced a new physical observable, "reaction". 

Define a complex number denoted by <a, b> as follows: 




V 




(ID 



where the volume V contains the source a. 



The reciprocity theorem states that 



<a, b> = <b , a> 



( 12 ) 



provided that all media are isotropic and that a and b can 
be contained in a finite volume. The scalar <^a, b> is 



a measure of the reaction (or coupling) between the sources 
a and b. Equation ( 11 ) is to be understood as a formula 
for the measure of reaction. It is obvious that Equation (9) 
may be restated as Equation ( 12 ) , where the inner product 
represents the complex integrals of Equation (9). Equation 
(12) is, therefore, the reaction form of reciprocity. 

These theorems will now be used to develop the "Piecewise 
Sinusoidal Reaction Matching Technique" (PSRMT) of J. H. 
Richmond. 
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IV. PIECEWISE -SINUSOIDAL REACTION MATCHING TECHNIQUE (PSRMT) 



A. INTRODUCTION TO PSRMT 

In this section, the method of moments and the reaction 
concept will be applied to the solution of the current 
distribution on linear antennas. This particular class 
of problem involves only antennas which consist of or can 
be represented by a collection of thin-wire segments. 

An axial filamentary test source will be developed. 

After specifying a suitable current distribution on both 
the original and the test antenna, the reaction concept 
will be applied. Reaction will be employed in order to 
develop a system of linear- equations whose solution will 
yield the complex coefficients that complete the description 
of the approximate current distribution specified. 

The reaction concept and its applications have been 
discussed by Rumsey [17] and Harrington [9]* The interpre- 
tation of PSRMT to follow is primarily that of Richards [13] 
and Leonard Tsai [21]. Most of the material presented 
here has been previously presented by Richmond [14, 15, 16 ] 
and Richards; it is presented here for completeness. 

B. THIN-WIRE PROBLEM 

Consider the thin-wire dipole shown in Figure 1. The 
dipole arms are wire segments of length L^ and , not 
necessarily equal, but shorter than ^/2 . The radius "a" 
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equivalent Surface Current 




FIG'JP.E 1. A thin -wire dipole wi in a fila- enter/ 
tost source on its e/is. 






of the wire is arbitrary, but in the thin-wire case always 
small in terms of A; also, a « and a « The 

dipole is driven at the junction of its arms with a magnetic 
ring current 

Let the dipole depicted in Figure 1 radiate in free 
space, setting up time-harmonic field intensities (EL., H^) 
everywhere. By the surface equivalence principle this 
antenna may be replaced by a set of surface currents 
(J d , M^) enclosing the primary source where 

J d = n X H d (13A) 

M d = E d X n (13B) 



at the surface, and n is the unit normal vector pointing 
outward from the surface. Let this surface exist at an 
infinitesimal distance outwardly away from the actual 
dipole, so that the dipole itself is replaced by surface 
currents (J d , M d ). Then it is known that exterior to this 
surface the two fields (E d , H d ) exist and there is a null 
field inside the surface. Now the original source can be 
removed without disturbing the fields, so that there is 
free space everywhere. 

A test source is placed at the axis of the actual 
dipole in Figure 1, in order to employ the reaction test. 
Alone, this filamentary test dipole radiates the fields 



2H 



(E^_, H^_) in free space. The Carson reciprocity theorem of 
Equation (9) may be applied between the test dipole and 
the unknown equivalent surface current on the original 
dipole to give: 



where S g is the surface of the original dipole as shown in 
Figure 1, t is the coordinate along the filamentary test 
dipole, and ( 1^ , M^) represents the currents of the test 
dipole. Due to the null field generated by the original 
source in the interior of S Q , where the coordinate along 
the filamentary test dipole is, the right-hand side of 
Equation (14) is identically zero. Therefore the reaction 
between the two sources is zero. 



hence the term "zero reaction condition". This zero 
reaction theorem was developed by Rumsey [173. Equation (14) 
written in light of the zero reaction condition is: 




e 



(14) 



<a, b> =• <b, a> - 0 




e 



e 



(15) 



For a typical antenna problem, the excitation consists 
of a voltage generator driving the dipole at the junction 



of its arms. Since the antenna is thin, a "delta gap" 
voltage generator is assumed. Here the impressed current 
densities are: 



J d = 0 (16A) 

M d = fS ( £ - a, l - l Q )V d , (16B) 



where is the strength of the voltage generator, "a" 
is the radius, l q is the value of the coordinate £ at the 
gap and 6 is the Dirac delta function. 

In this case M d is a ring of magnetic current at the 
terminals, and the integral on the left-hand side of 
Equation (15) reduces to a line integral. Applying Ampere's 
law to perform the integration around the gap: 



If 



r2i r 



M d • H t ds e = -a V d 



J H t (a,£ 0 ,*5)d^ 



= ' V d X t (£ 0 } * 



(17) 



Combining Equations (15) and (17) yields 

Jf 1 d ' E t as e ’ - V d h ( V (18) 

b e 

where 1^ (£q) is the complex terminal current on the test 
dipole . 
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For the evaluation of Equation ( 1 8 ) allow the thin-wire 



approximation (a « "M to hold, thus the distribution of 
current density may be assumed to be uniform around the 
cylinder. This assumption may be expressed mathematically 
as simply, = I^/2ira where "a" is the radius and 1^ 
varies only along the length of the dipole. Thus Equation 
( 1 8 ) becomes 



Equation (19) is the result of replacing the surface 
current density with an equivalent filamentary current 
of strength X ^ = 2-rraJ^ which is located on the previous 



this line on the dipole surface and E^. represents the field 
of the test dipole on this line. It should be noted that 
if the test dipole were not axially located, then this 
simplification would not be judicious. 

Equation (19) forms the foundation upon which the two 
segment reaction matching technique is based. Once 
is specified and the current distribution 1^ is assumed, 
the coefficients for the current X^ may be readily determined. 
For the remainder of this paper the discussion will involve 
current distributions on both the test dipole and on the 
original dipole that are "piecewise-sinusoidal" in form. 

This will lead to considerable simplification because 



V 



d 




(19) 



e 



surface where resides. The integration d£. g is along 
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a.) closed form expressions for the E^_ term can be found, 
and b. ) closed form expressions for \\ 1 ' E t W e also 
exist for coplanar dipoles. The basis for these claims 
will be presented after a generalization is made for N 
segment reaction matching. 

An essential feature of any acceptable test source 
(J^., M^) must be that the fields it radiates be known and, 
if possible, that they be known in closed form. One such 
source is a filamentary electrical source with time-harmonic 
current, as shown in Figure 2. 

Let the source be located on the z axis with endpoints 
z^ and z^, and the current distribution be of the form: 



I(z) = A e 



jkz 



B e 



-jkz 



( 20 ) 



where A and B are complex constants, k = 2 tt/x, and the 

1 wt 

time dependence e^ is understood and suppressed. 



C. ELECTRIC FIELDS OF A PIECEWISE-SINUSOIDAL FILAMENTARY 
SOURCE 

As expression for the z component of the free-space 
electric field of the filamentary source of Figure 2 is 
given by Schelkunoff and Friis [19] as 



1 

i|-TQ U) £ 



I i 



,-jkR. 



1. T I ® 

R„ 2 R 



-JkR, 



Here I. = I(z n ), I* = I'(z n ), I 0 = I(z 5 ) 



>V 

>z 



( 21 ) 



1' 2 = I’(z 2 ), 
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Figure 2 







• origin 



Filamentary electric line source on the z axis 
sets up a field (E f ,E ) at an observation 
point U,z). 
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and V is the electric potential of the end charges; £ 
is the permittivity of the medium. Throughout this paper 
the prime denotes taking the derivative with respect to z. 

The electric field is composed of both a normal 
component and a tangential component E as depicted 
in Figure 2. 

Note that in performing the dot product "J^ • E^" 
of Equation ( 1 8 ) for linear antennas, only the z component 
of the electric field makes a contribution. That is the 
justification for ignoring the E^ component when considering 
only linear segment configurations. 

The term ^V/^z will be dropped immediately since it 
vanishes for adjacent segments and the contributions from 
the gradient of the potential V due to end charges for the 
end segments may be assumed negligible [13]. Note that 
although this expression for E is relatively simple, this 
closed form expression is rigorous even in the near-zone. 

The complex coefficients A and B in the expression 
for I(z) (see Equation 20) can be solved in terms of the 
end currents 1^ and I 2 * The resulting equation for I(z) 
is 

I(z) = [1^ sin k(z£ - z) + sin k(z - z^)] /sin kd 

( 22 ) 



where d as indicated in Figure 3 is the length of the 
segment . 
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FIGURE 3 
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Piecewise-sinusoidal current distribution 
on an antenna segment. 
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Figure 3 illustrates the current distribution expressed 
in Equation (22). It is clear that from Equation (22) 

I£ and I), may also be expressed in terms of 1^ and I 2< 

Since 

I’(z) = C-kl^ cos k(z 2 - z) + kl 2 cos k(z - z^)] /sinkd 

(23) 



it follows that 

1^ = k(I 2 - I 1 cos kd) / sin kd (2*JA) 

IJ = k(I 2 cos kd - I ) / sin kd . (2*JB) 

Note that from Equations (21) and (2*0 the electric field 
intensity of the sinusoidal electric line source may be 
expressed entirely in terms of the geometry of the observation 
point and the end currents 1^ and I 2> 

D. EXPANSION OF THE CURRENT DISTRIBUTION 

The forms of the current distributions on both of the 
filamentary dipoles is piecewise-sinusoidal . The current 
on each segment of each dipole can be described by the 
sinusoidal functions given in Equation (22). Because the 
components of the electric field intensity of such an 
element are known, both terms within the integral of Equation 
(19) exist as closed form expressions. 



32 



The following piecewise-sinusoidal distribution 



has been employed by Richmond [153: 



J 





sin k {l-l .. ) 
n-1 



s 



sin kd 




sln kU tvfl - l }_ (25) 



sin kd 



or 



J 



y 1 o 

2tta n n 



1 



( 26 ) 



s 



with I = I(£ ) and d the length of the segment. P (£) 



is a unit pulse, nonzero only on segment n; T is a unit 
vector pointing in the direction of l from one endpoint of 
segment n to the other. The vector functions are 
known as basis functions, expansion functions or dipole 
modes. The mode currents go to zero at the endpoints 
of the dipoles, as depicted in Figure 4. 

It is advantageous to have a symmetric impedance matrix 
Z mn (to be defined later). Furthermore, the test sources 
should be chosen to yield a well-conditioned system of 
linear equations for the subsequent inversion process. 

For these reasons and to obtain a closed form expression 
for the integral of Equation (19), the expansion function 
for the test source G„ will be of the same form as the 



basis function for the original dipole's current distribution 



m 



G. 



n* 
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l n-Z Vl £ n+l 



FIGURE 4. Basis functions (G n ) representing the 
assumed current distributions on a 
linear filamentary antenna. 
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With the piecewise-sinusoidal current distribution 



employed for each segment of the antenna, then from 
Equation (22) the current on any given segment can be 
considered the sum of two sinusoidal currents, each of 
which vanishes at one or the other end of the segment. 

This may be visualized with the aid of 'Figure *1. One can 
see how the equivalent surface distribution is treated as 
a collection of (N-l) overlapping tubular dipole current 
densities on an antenna consisting of N segments. The 
nth dipole has terminal current I at the nth junction with 
the current vanishing at the far ends of the two segments 
forming the nth junction. 

The solution for the unknown complex terminal currents 
( I n> proceeds by introducing a set of (N-l) filamentary 
dipoles each having piecewise-sinusoidal currents and 
being located along the axis of the original antenna. No 
generality is lost by assuming each test dipole has unit 
terminal current. Reciprocity is enforced successively 
between each test dipole and all the tubular surface 
dipoles in turn. The "zero reaction condition" yields 
the following equation for the mth test dipole and the 
nth surface-current dipole (refer to Equation (15)): 



The left-hand side of Equation (27) is nonzero only 
when n = m. This is the result of the tangential component 




(27) 



n 



n 
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of the electric field on the surface of the original source 
being zero except at the delta gap. which equals 

E x n is therefore also zero unless the mth test dipole 
is co-located with the nth original dipole. Expressed 
mathematically : 






d l 

n 



mn 



-v i £ 

n m mn 

( 28 ) 



The right-hand side of Equation (27) can be simplified 
by carrying out the circumferential integration of • E . 
Again this is equivalent to replacing the tubular distri- 
bution with filamentary segments for calculating reaction: 



s 



n 



m 



ds 



n 



f T U ) 

n' 



E d l 
m n 



(29) 



n 



n 



Equating Equations (28) and (29) one obtains 



-V I i = fl (J ) • E d^ 
nmmn PJ n n m n 

''n 

Assuming unit amplitude for the test current and noting 
that upon co-locating m = n. 



V <f = - fl (t )-E dT (31 

m mn J n n m n 

v n 

The integration is carried out on^ n> the coordinate along 
the arms of the nth surface dipole filaments. The 
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symbol & mn is the standard Kronecker delta, defined by 
the relations: 



£ =0 m ? n 

ran 

7 * 0 m = n 



( 32 ) 



Equation (31) forms the basis of the multi-segment reaction 
matching technique. 



E. MATRIX EQUATION FORMULATION 

Substitute the assumed form of I n (+- ) (refer to I'igure 
*0 namely 



I (r ) 

n n 



= I 



n 



sin k(l - l , ) 

£ P U) 7 — r t ?- ” - 1 ~ 

n n sin kd 



sin k (l ,, 

+ £ P u) nil. 

n+l n+l''*"’ sin kd 



- 1 ) 



(33) 



into Equation (31) and define a mutual impedance term as 



mn 



■ 



~n-l 



“n sin k (l - l , ) 

n -~- T • E dZ 

sin kd n m 



J l n+l si 



nk(£ , - l) ^ 

— 'b • E d£. (3*0 

n+l m 



n 



sin kd 



The dummy variable l is the coordinate along the particular 
segment, the other variables are defined as in Equation (25). 
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E m as before, is the free space electric field intensity 
radiated by the mth test dipole with piecewise-sinusoidal 
current, specifically for the mth test dipole composed of 
a lower segment ML and an upper segment MU, the z component 



of the electric field is 
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where as 


depicted in Figure 5, 
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Figure 5 is merely an extension of the single segment 
with piecewise-sinusoidal currents of Figure 2 to a dipole 
composed of two such segments. 



Equation (31) may now be written 



V 

m 





(36) 



where I is the complex terminal current at the nth junc- 
tion. The calculation of the reaction between the mth 
test dipole and the entire equivalent surface current 
density, consisting of N-l overlapping surface dipoles, 
yields the following relation: 
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origin 



FIGURE 5. 



(*.z) 




f'h 

The m Lr test dipole, composed of a lower segment 
ML and an upper segment MU 
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V 

m 



(37) 



N-l 

Z Z 

n=l 



mn n 



Equation (37) represents one linear equation out of a 
system of (N - 1) equations that are obtained by successively 
enforcing reciprocity (or calculating reaction) between the 
(N - 1) surface dipoles and each of the (N - 1) axial test 
dipoles in turn. The test is performed in turn between 
each of the axial dipoles and the tubular dipoles yielding 
a system of linear equations that can be described by the 
matrix relation 



[Z ] 
mn 






tv ] 

m 



(38) 



which is of the same form as Equation (37), except that 
Z mn is now a square (N-l) X (N-l) matrix and V m is a 
column vector. This set of equations can be solved for 
the unknown junction current coefficients -^1^ after a 
known voltage source is placed at the actual feed junctions 
and the voltage of the other junctions are set to zero. 

The coefficients I n of Equation (38) can be solved for 
by several techniques. Crout published a paper [7] in 
which he discusses a method for evaluating determinants 
and solving systems of linear equations with complex 
coefficients. In the case of well conditioned matrices 
the impedance matrix can be merely inverted and multiplied 
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by the voltage matrix to yield the terminal current 
coefficients : 






[z ] 

mn 



-1 



[V ] 
m 



Equation (39) may be written simply as - 

-1 

[I] = [Z] [V] 



( 39 ) 



(40A) 



or equivalently: 

[I] = [Y] [V] , (40B) 

where 

-1 

[Y] = [Z] 

The complex current coefficients 1^ that are determined 
from Equation (38) are used to establish the current 
distribution found from Equation (26). Once the current 
distribution is known, parameters of interest such as field 
patterns, input admittances, etc. , can be calculated by 
numerically evaluating the appropriate formulas. 

A wire antenna is obtained when the wire is excited 
by a voltage source at one or more points along its length. 
For an antenna excited at the nth junction, the applied 
voltage matrix of Equation (40) is 



0 




V 



n 



(4l) 



0 



that Is, all elements of the voltage matrix are zero except 
the nth, which is equal to the source voltage. The current 
distribution is established by Equation (40) which for 
the [V ] of Equation (4l) becomes 



Note that the nth column of the admittance matrix gives 
the current distribution for a unit voltage source applied 
to the nth junction. Inversion of the impedance matrix 
therefore gives simultaneously the current distributions 
when the antenna is excited at any arbitrary junction 
along its length. For the case where only one junction 
is excited with a voltage , the input admittance is 
given directly by the feed junction current 1^ as 




(42) 
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P. SOLUTION TECHNIQUES 



As mentioned previously, integrals of the form of 
Equation (3^) can be obtained in closed form for coplanar 
segments. In the case of non-coplanar segments, the 
integrals must be calculated by numerical methods, for 
example, Simpson's rule or Gaussian quadrature. 

Upon a careful examination of Equations (3^0 and (35) 
it should be ovbious that the integral of Equation (35) 
may be expressed in closed form with the relationship [16] 
that follows: 
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(44) 



where 



Ci ( x) = V + Ln(x) + J dt (45A) 
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Si (x) 



t 



dt 
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(45B) 



Utilizing the relationship in Equation (44), Richmond 
and Geary published a paper [ 1 6 ] that developed a closed 
form expression for the mutual impedance between coplanar 
skew dipoles. However, for this antenna problem, all of 
the reactions will be performed between parallel dipoles. 

A more useful closed form expression for the problem 
consisting of entirely parallel dipoles was presented by 
H. E. King [10]. The complete expression with its deriva- 
tion may be found in Appendix A and B. Even though these 
closed form expressions are not particularly simple, they 
do involve only sine and cosine integrals which have 
polynomial algorithms; hence Z in this form is more 



efficient than "brute-force" integration. 

Leonard Tsai [21] formally portrayed PSRMT as a 
particular application of Galerkin's method. Consider 
again the operator equation 



where f, the given operator, is specifically the unknown 
current distribution on a filamentary dipole and the 
operator L yields the fields g of a distance "a" away 
(on the original dipole). Expanding f into a linear 
combination of basis functions: 



L(f) = g 



(46) 



f 



n 




(47) 



and enforcing Galerkin's method directly: 




( 48 ) 



Here the f would describe the weighting functions 
corresponding to the "test dipoles" which is using the 
same expansion functions as the "original dipoles", both 
subdomain piecewise-sinusoidal pulses. The inner product 
of the left-hand side of Equation (48) corresponds directly 
to the integral of Equation (31): 



used in the reaction matching technique. Note also that 

<g, f> = V which is nonzero only at voltage sources. 
' mm 

This follows from the "zero reaction condition". Thus 
the procedures of PSRMT are those employed by a direct 
application of Galerkin's method. 




(49) 
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V. DISCUSSION OF NUMERICAL TECHNIQUES AND COMPUTER PROGRAM 



Evaluation of the integrals in Equation ( 3*0 may be 
accomplished by any one of several numerical integration 
techniques such as Romberg, Weddle, Gaussian Quadrature, 
etc. , as well as by the closed form expression developed 
in Appendix A. For most integrals, no representation in 
terms of elementary functions is possible, and approximation 
becomes necessary. Most numerical techniques are successfully 
based upon some type of polynomial approximation: trapezoi- 

dal rule, Simpson's rule, Romberg's method, etc.. Unfortu- 
nately, these techniques have the usual error sources such 
as truncation, but perhaps more interestingly they are 
harassed with the question of convergence. The question 
is whether as continually higher degree polynomials are 
used or as continually smaller intervals between data 
points are used, a sequence of approximations is produced 
for which the limit of the truncation error is zero. As 
may be guessed, the smaller the interval h is made, the 
more computation time required and the greater the rounding 
off error becomes. In other words, accuracy is limited 
for any given method by algorithm errors that ultimately 
obscure convergence. 

It is found in practice that decreasing the interval 
between sampling points below a certain level leads to 
larger errors rather than smaller. This necessitates 
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experimentation of the various parameters available in 
the technique in order to establish a trend (if one exists) 
that gives acceptable results. Unfortunately this requires 
an assumed correct solution in order to evaluate the 
performance of the integration on a supposedly unknown 
integral. 

Table 1 illustrates how the value for input admittance 
of a half-wave dipole varies as a function of the number 
of antenna segments and the order of the Gaussian Quadrature 
used in evaluating the integral of Equation (3^). Table 2 
is a similar portrayal of a full-wave antenna. 

The norm for Table 1 and 2 was chosen to be the value of 

input admittance given by King [11]. Table 1 shows no 

convergence tendency whatsoever, while Table 2 indicates 

that a lower degree polynomial may be used in evaluating 

Z if the antenna is divided into more segments. The 
mn 

apparent discrepancy between Table 1 and 2 is probably 
explained by the difference in the actual current distri- 
butions on the two antennas. It is clear that the problem 
is multi-faceted. 

This convergence problem is not exclusive to the 
Gaussian Quadrature technique. In fact, all of the 
numerical techniques have some troublesome parameter analysis 
involved in order to effectively use them. But a far more 
important and exciting reason to avoid their use is that 
of efficiency. The closed form expression involves only 
the evaluation of the sine and cosine integrals which have 
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GAUSSIAN QUADRATURE 




TABLE 1. Convergence analysis for the input admittance 
of a half-wave dipole. The numbers respresent 
percent variation from the King-Mi ddleton standard 
of Y= 10.8 mmho's. (Percentages expressed to the 
nearest .5%) 
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GAUSSIAN QUADRATURE 




TABLE 2. Convergence analysis for the input admittance of a 
full -wave dipole. The numbers represent percent 
variation from the King-Mi ddl eton standard of 
Y = -SO mmho's. (Percentages less than 20 expressed 
to the nearest .5%) 
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simple algorithms. Unfortunately the closed form solution 
is valid in general for skewed coplanar dipoles only, 
therefore, numerical integration would be required -for 
non-coplanar elements. The same antenna problem (if 
coplanar) may be solved by the closed form technique in 
a fraction of the time that numerical methods require. 

The computer program listed in Appendix B calls upon 
subroutine CZMN2 for the evaluation of the impedance 
matrix. It is here that the method of solution may be 
chosen. If the closed-form route is decided upon, 

CZMN1 (which is also listed) is used to evaluate Z . 



v 
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VI. ANALYSIS OF RESULTS 



A. CONVERGENCE CHECK 

Figures 6 and 7 illustrate the accuracy and efficiency 
of PSRMT. Leonard Tsai [21] used similar figures (the 
solution labeled "Richmond's Reaction") to show the compara- 
tive rate of convergence of PSRMT to that of point matching 
with triangular pulses. Note that reaction matching 
required only six segments while for the same accuracy 
point matching required five times as many segments. 
Superimposed on Tsai's comparison are those points generated 
by the closed form expressions of Equation (72) of Appendix 

B. 

Figures 8 and 9 illustrate the efficiency of PSRMT. 

In Figure 8, as more and more flat pulses are employed 
in Galerkin's method the more the solution corresponds to 
that of PSRMT with only six segments. Figure 9 shows how 
rapidly PSRMT relates to the moment method employing 
piecewise-sinusoidal basis functions. 

Figure 10 exemplifies the convergence of PSRMT itself. 
The comparison is drawn between PSRMT and Hallen's equation 
employing piecewise-sinusoidal basis functions. 16 segments 
do an admirable job in relating to this half-wave dipole. 
This is sufficiently true under certain conditions (as 
will be shown) to warrant a rule of thumb: use 32 segments 

per wavelength. 
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Figure 11 is a plot of conductance versus susceptance 
for thin cylindrical dipoles with radius a/\ = .001588. 

The locus spirals in as the length of the dipole increases. 
Note that for small lengths the 8 segmented solution devi- 
ates from the accepted King-Middleton value of admittance. 
Many of the plots to follow will use King-Middleton values 
for comparison, which were taken from tables published by 
R. V/. P. King [11]. In the evaluation of the admittance, 
the King-Middleton procedure was again used (two iterations 
and an empirical normalization were also employed) . Figure 
11 indicates the advantage of taking 32 segments per 
wavelength. 

Figure 12 is the same type as Figure 11 but the radius 
has been increased to .006351. It is evident that for short 
dipoles with large radii, the length to radius ratio for 
each segment exceeds the thin-wire approximation. That is 
why that as the length of the antenna is increased in Figure 
12, the solution becomes more convergent. 

Figures 13 and 14 portray the convergence trends of the 
Weddle numerical integration technique on a half-wave dipole. 
As indicated in Figure 13, the solution varies only slightly 
with a change in the number of integration points used in 
the integration. However, the input admittance approaches 
the King-Middleton solution most nearly when the number 
of integration points equals 29. Figure 14 shows the change 
in current distribution as the number of antenna segments 
change (with the number of integration points equal to 29). 
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Again the input admittance approaches the King-Middleton 
solution as the number of segments increase. These 
particular trends are valid only for this half-wave dipole 
of radius a/x = .001588. Other geometries must be considered 

on a one to one basis due to the behavior of Weddle's 
technique on the integral of Equation (3*0. 



53 



R 

(ohms) 




FIGURE 6. Input resistance versus h/* for a cylindrical wire 
antenna with radius a and total length L = 2h. 
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FIGURE 7. Input reactance versus h/-^ for a cylindrical wire 
antenna with radius a and total length L = 2h. 
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Gal erki n Method with Flat Pulses 



6 pulses 
13 pulses 
20 pulses 




FIGURE 8. Current distribution on a half-wave dipole with 
radius a/*. = .001 588. 

(Galerkin Method with flat pulses) 
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General Moment Method 



I (ma) 



if 



6 basis functions 

*13 basis functions 

20 basis functions 



30 basis functions 

PSRMT - Closed Form 




FIGURE 9. Current distribution on a half-wave dipole with 
radius h/>. = .001588. 

(General moment method — piecewise-sinusoidal pulses) 
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Hallen's Equation 

Symbol Divisions Input Impedance 




FIGURE 10. Current distribution on a half-wave dipole with 
radius a A = .001. 

(Hallen's Equation - piecewise - sinusoidal - basis functions) 
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King - Middleton 
PSRMT - Closed Form 



A 8 segments 

o 32L/% segments 



B 

(mmho's) 




FIGURE 11. Admittance curve for thin-cylindrical dipoles with 
radius a/% = .001588. 
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King - Middleton 



PSRMT - Closed Form A 8 segments 



+ 16 segments 
B 24 segments 




FIGURE 12. Admittance curve for thin cylindrical dipoles 
with radius a /x = .00635. 
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(numerical integration) 5 
Weddle's Method 17 

6 segments) 29 
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FIGURE 13. Current distribution on a half-wave dipole with 
radius a/v = .001588 (convergence check on the 
number of integration points used in 
Weddle's technique) 
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Symbol Y 




FIGURE 14. Current distribution on a half-wave dipole with radius 
a/x = .001588. 

(convergence check on the number of segments used in 
Weddle's technique) 
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B. ADMITTANCE CORRELATION 

Figures 15 through 18 attempt to typify the admittance 
correlation between the King-Middleton solution and various 
PSRMT solution techniques. Figure 15 and 16 depict the 
solution by Weddle's method and the Gaussian Quadrature 
method, respectively. Figures 17 and 18 can be compared 
with Figure 11 to show the effect of changing the dipole's 
radius. The effect is negligible as long as the length 
to radius ratio of each segment is maintained large. 
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FIGURE 15. 



PSRMT - Numerical Integration ^ 
(Weddle's Method 
Number of Integration 
points equals 29) 




Admittance correlation for thin cylindrical dipoles 
with radius a/->. = .001588. (32 L/;*. segments) 
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Ki ng-Middl eton 



c 



PSRMT - Numerical a 

Integration 
(Gaussian quadrature 




FIGURE 16. 



Admittance correlation for thin cylindrical dipoles 
with radius a/x = .001588. (32 L/>. segments) 
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FIGURE 17. Admittance correlation 
radius a/> = .003175. 



for thin cylindrical 
(32 L/-v segments) 



dipoles with 



King-Mi ddleton * 

PSRMT - Closed Form a 

B 

(mmho ’s) 




FIGURE 19. Admittance correlation for thin cylindrical dipoles 
with radius a/j. = .004763. (32 L/* segments) 
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C. CURRENT DISTRIBUTIONS 



Figures 19 through 25 present current distributions on 
dipole antennas of various lengths and radii. On each set 
of curves is indicated the expansion set used to represent 
the unknown current. The number of segments employed in 
the PSRMT technique was based upon the rule of thumb 
32 LA . 

In Figures 24 and 25, the power series [2] used as a 
comparison may be expressed as 

I(Z) = Jfi Zn " 1 1 

where the C 1 s are complex constants to be determined. 
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King-Mi ddleton 



PSRMT - Closed Form A 
(8 segments) 




FIGURE 19. Current distribution on a quarter-wave dipole with 
radius a/% = .001589. 



69 




FIGURE 20. Current distribution on a quarter-wave dipole with 
radius a/% = .007022. 
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King-Mi ddleton 




FIGURE 21. Current Distribution on a 3/4-wavelength 
dipole with radius a/7- = .001588. 
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King-Mi ddleton 




FIGURE 22. Current distribution on a full -wave dipole 
of radius a/x = .001588. 
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FIGURE 23. Current distribution on a 3/2 > dipole of radius 
a/>. = .001588. 
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King-Mi ddleton 
Butler-Power series 
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F 1 GU RE 24. Current distribution on a half-wave dipole with 
radi us a/>. = .001 588. 



Butler-Power series 
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(16 segments) 

- Gaussian Quadrature 
(4 segments) 
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FIGURE 25. Current distribution on a half-wave dipole with 
radius a/}. = .005. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A. EVALUATION OP POTENTIAL INTEGRALS 



Integrals of the form 




f(z') [- 





e 



] dz 1 



( 50 ) 
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1 



R 



2 



where 
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1 




R 
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~\J~( z+z ’ ) 2 + 



are extremely difficult to evaluate for "small a" because 
of the sharp peak in the integrand at z' = z. Although 
PSRMT has been formulated on the thin-wire assumption 
(a <<^); for very small radius values the integrand in 
Equation (3*0 becomes ill-behaved due to the form of the 
electric field E m (refer to Equation (35)). 

This phenomena was observed in the calculation of the 
input impedance of Figures 6 and 7- In fact, for h/w <( .3 
the value of the radius (h/a = 2000) became small enough 
to thwart computation entirely. 
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B. APPLICATION OF THE GENERAL PSRMT TECHNIQUE 



Although this paper dealt in particular with parallel 
dipoles only, PSRMT is useful in solving a wide variety 
of thin-wire antenna problems. With the aid of the concepts 
of reciprocity, reaction and surface equivalence, this 
technique reduces an integral equation to a set of simul- 
taneous linear equations. This set of equations is readily 
solved on a computer to obtain the input impedance and 
current distribution on the antenna. 

Richmond and Geary [ 1 6 ] presented a new form for the 
rigorous near-zone field of the linear dipole with sinusoidal 
current distribution. This electric field expression was 
employed to derive a closed form expression for the mutual 
impedance between coplanar-skew dipoles. By employing 
their expression for mutual impedance, complex geometries 
and surfaces may be represented by wire grids and reaction 
enforced to yield solutions. 

C. RECOMMENDATIONS 

The theory developed in this paper is by no means 
restrictive to the specific problem investigated. The 
data comparisons were made by employing PSRMT to equally 
segmented, thin-wire linear, parallel dipoles. Extensions 
to this problem are obvious. First, the segments need not 
be of equal length. By merely altering the computer program, 
the specific segment lengths can be applied to the appro- 
priate formulas to generate a more general solution. Second, 
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the dipoles need not be parallel. By employing Richmond 
and Geary's closed form expression, the mutual impedance 
between any two coplanar elements may be determined. 

But now the reaction between the filamentary test dipole 
and the equivalent filamentary tubular surface dipole must 
be examined in more detail. In the case of parallel 
elements, the reaction applied between the test dipole and 
any filamentary surface dipole is the same. But due to 
the geometry of the problem involving two non-parallel 
segments, the mutual impedance depends upon the angular 
position of the filamentary surface dipole with respect 
to the axis. The question of how does one best react these 
segments is yet to be answered. E. G. Neely will present 
a Master's Degree Thesis to the Electrical Engineering 
Department at the U. S. Naval Postgraduate School (in 
September 1973 ), in which he deals with arbitrarily bent 
monopoles above a ground plane. He is approaching this 
problem by averaging the reactions taken around the 
circumference of the tubular dipole successively. Other 
approaches might be investigated. 
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APPENDIX A 



TRANSFORMATION OF EXPONENTIAL INTEGRALS 
INTO CLOSED FORM EXPRESSIONS 

The z component of the electric field from a dipole 
with a time-harmonic current distribution on each segment 
is expressed by Equation (35)- For the mth test dipole, 
the dipole mode current distribution is shown in Figure 26. 
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FIGURE 26. Typical mode current distribution for a 
filamentary test dipole. 
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